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Abstract 

A characteristic action AS" is defined whose magnitude determines some properties of the ex- 
pectation value of a general quantum displacement operator. These properties are related to the 
capability of a given environmental 'monitoring' system to induce decoherence in quantum systems 
coupled to it. We show that the scale for effective decoherence is given by AS" ~ h. We relate this 
characteristic action with a complementary quantity, AZ, and analyse their connection with the 
main features of the pattern of structures developed by the environmental state in different phase 
space representations. The relevance of the AS-action scale is illustrated using both a model quan- 
tum system solved numerically and a set of model quantum systems for which analytical expressions 
for the time-averaged expectation value of the displacement operator are obtained explicitly. 

PACS numbers: 03.65.Yz, 03.65.Ta 
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I. INTRODUCTION 



The superposition principle and the interference terms that generates are the key com- 
ponents of the quantum formahsm, and responsible for the main differences between the 
quantum and classical world. The boundary between these two worlds and the mechanisms 
that prevent the interference terms from being apparent in the classical realm have been 
the subjects of many theoretical and experimental studies since the very beginning of the 
"quantum era". Significant advances in the analysis and experimentation on the interac- 
tion between mesoscopic and microscopic systems are pushing the boundary between the two 
worlds. An example is the study of measurement processes where the 'monitoring' apparatus 
is represented by a system with an increasingly larger number of degrees of freedom (more 
classical) and the analysis of the associated disappearance of the non-diagonal terms of the 
density operator of the microscopic system in some preferred matrix representation 
The study of the effectiveness of a given system that plays the role of an environment or of 
a measurement apparatus to induce decoherence in another system is of fundamental and 
practical interest. For instance, the advances in the fields of quantum communication and 
quantum computation depend crucially on our ability to manipulate entanglement and 
to control the capability of the environment or measurement devices to induce decoherence 
in our qubit (pointer) system [J, 0|. 

Many actual interactions between a two-level system S, spanned by the pointer states |-|-) 
and |— ), and a system S playing the role of the environment (for instance as a 'monitoring' 
apparatus), can be described by means of a coupling Hamiltonian of von Neumann's form 
j^. In particular, we will use a generic term V^g = (|+) (+| — | — ) (— |) (Cq ■ q + Cp ■ p), 
where q = (gi, . . . , g/) and p = {pi, . . . ,pf) are position and momentum operators for an 
environmental system with / degrees of freedom {[qj,Pj] = ih, j = 1, . . . , f). The coefficients 
Cq = {c'q\ . . . , c[/'') and Cp = {c^\ . . . , c/^) characterise the strength of the coupling. The 
reduced density operator describing the state of the system S after its coupling with the 
environment during a time interval 6t is given by 

PS = i«n+)(+i + i/?n-)(-i +(«/?* (^-iv^+)i+)(-i+H.c.) , (1) 

where \iIj±) = D{^Cp6t, =FCq5t) lip), -D(5q, 6p) = exp{z(p ■ 6q + q- 6p)/h}, 5q and Sp are 
displacement vectors in /-dimensional spaces, and H.c. denotes the Hermitian conjugate of 
the preceding term in the equation. The states of the environmental and two-level systems 



immediately prior to the interaction are \ip) and |x) = « |+) + /^|~) respectively. We 
have assumed that the coupling strength is large enough so that the evolution induced by 
each system Hamiltonian [Hs and Hs) can be neglected during the interaction time 5t j^. 
Despite the simplicity of the model considered, it contains the basic elements relevant to 
our discussion. 

Eq. (P) relates the value of the non-diagonal term of the reduced density matrix of S in 
the preferred basis {|+), |— )} to the mean value of a displacement operator over the state 
I?/') of system £ since {ip_\ip+) = D{—2Cp6t, — 2Cq5t) Therefore the capability of S 
to induce decoherence in S through the coupling term Vss is characterised by 

where ip{q) = (q|^), 5q = —2Cp6t, 6p = — 2Cq(5t, and d^q (d^p) is the /-dimensional 
differential element of volume in positions (momenta). All integrals in this paper run over 
the entire available volume. Complete decoherence is reached whenever the two states 
and lip-) are orthogonal to each other; in other words, when = 0. At this point, it 
is important to characterise the scale for which displacements (5q, 6p) in phase space will 
produce a significant decay of this expectation value of D. The main subject of our interest 
is to find an action scale associated to the effectiveness of system S to induce decoherence 
in system S, and to describe its dependence with the particular environmental state. 

space distribution associated to the state lip) p, 

W4ci,p) = I rfVe^'^'-P/'^(q-q72)V^*(q + q72). (3) 

In particular, Moyal's formula Ql 

\C46q,5p)f = {2nhy J d^qdfpW4c^,p)W4ci+6c^,p + 5p) (4) 

was used to analyse the behaviour of the overlap |C^P with 5q and 6p. The choice of 
the Wigner phase space distribution was motivated by this simple expression for the scalar 
product between | ■?/'+) and IV'-)- In Ref- Zurek showed that for a given time- dependent 
quantum chaotic system in one dimension (/ = 1) confined to a phase space volume charac- 
terised by the classical action A, the Wigner distribution associated to the state develops in 
time a spotty random structure on the scale h'^/A. Using Eq. (jlj he argued that |C^P ~ 



for phase space displacements on the scale of the smallest structure of the Wigner distribu- 
tion W^{q,p). The basis for this result are: (a) Displacements characterised by 6q6p ~ h'^/A 
produce a significant decrease on the value of the integral in Eq. (jH) due to the destructive 
interference between W^{q,p) and W^{q + 6q,p + 6p), and (b) the random distribution of the 
patches in the structure appearing in the Wigner function associated to such system states 
prevents the presence of recurrences in the value of the overlap. 

Jordan and Srednicki extended the analysis in Ref. ^ to systems with an arbitrary 
number of degrees of freedom by using 

This equation establishes a relation between the small-scale (large-scale) structure of in 
the variables (q, p) and the large-scale (small-scale) structure of in the variables (5p, 6q). 
Analysing a two-dimensional billiard and a gas of hard spheres in a three-dimensional box 
(assuming the Berry- Voros conjecture in both cases) they concluded that for systems with 
a small number of degrees of freedom, displacements 6qi ~ Lj and 6pi ~ Pi are needed to 
avoid oscillations in the overlap, where Lj and Pj are typical classical values of the position 
qi and momentum pi respectively {i = 1,...,/). This means that displacements of the 
order of the size of the state support are needed to guarantee orthogonahty in the general 
case. However, for systems with a large number of degrees of freedom they found that the 
conclusions in Ref. ^ remain valid, supporting the idea that a larger number of degrees of 
freedom increases the effectiveness in causing decoherence. Some care must be taken when 
relating the results in Refs.S] and jsf since in principle the Berry- Voros conjecture is not 
valid for the system analysed by Zurek in Ref. P| and the dependence of the overlap with 
the displacement could have qualitatively different features. 

In this work we characterise the behaviour of using a quantity AS*, with units of 
action, associated to the displacement ((5q, 6p). A formal series expansion of D will allow us 
to identify the scale in the action AS'(5q, 6p) for which the overlap decreases significantly 
for any quantum system, irrespective of the number of degrees of freedom. This scale is 
manifested in the size of the structures present in the distribution associated to the state in 
some phase space representations, but they do not necessarily coincide. 

The paper is organised as follows. In Sec. II we define the characteristic action AS* and 
determine the scale relevant for the decay of the overlap. In Sec. Ill we establish the relation 
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between and the structure of the distribution associated to the state in an arbitrary phase 
space representation. The next two sections are devoted to studying in detail the dependence 
of on AS for states of particular quantum systems. Sec. IV considers a system with 
a time-dependent Hamiltonian whose classical counterpart exhibits chaos. In Sec. V we 
analyse the case of non-linear systems with a confining potential and discrete spectrum. In 
this case the main features of can be obtained from time average properties of the state 
evolution. We will focus on quantum systems with time-independent Hamiltonian for which 
analytical models are worked out by using the Berry- Voros conjecture jsl. Finally in Sec. 
VI the main results of this work are discussed. 



II. CHARACTERISTIC ACTION SCALES FOR THE DECAY OF C^,(5q,(^p) 

A displacement operator Z)(5q, 5p) acting on the state of an /-dimensional quantum 
system that describes an environment or a 'monitoring' apparatus, can be written as 

i)(5q, 5p) = e^^P-^'i+^-^P)/^ = e^^^'^'i'^P)/'' . (6) 

The main features of |C^P are therefore related to the fluctuation properties of the operator 
S'((5q, 5p), since the expectation value of D equals the characteristic function of S (see Eq. 
0). 

A formal expansion of D in terms of S gives 

and for the overlap, 

|C^(5q,5p)P = • (8) 

We denote by 5" general products of n components of the vectors 5q and 5p, and by 
terms of the form YYk^i{Ok)^^ where Ok is the product of operators q and p, with the 
condition ^'^^ Qk = The characteristic action 

AS{Sq, Sp) ^ ^J{S^)^ - {S)l (9) 

controls the decay of the overlap for sufficiently small values of Sq and Sp. Eq. (8) suggests 
that displacements {6q, 6p) for which AS* is small compared to h do not lead generally to 
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an important decay of |C^p. In other words, displacements leading to AS* of the order or 
larger than h are needed for the states and to be orthogonal. Therefore AS* ~ h 
establishes the scale for the action involved in displacements of the environmental state that 
could induce significant decoherence in system S. For the case of Gaussian fluctuations of 
the operator S, the only relevant fluctuation is AS*. In a more general situation higher order 
fluctuations may play a role in the particular features of the decay of |C^P, nonetheless the 
AS-action scale is generally expected to be a good measure for the decoherence process. 
The rest of the paper will provide additional arguments for this interpretation of the scale 
associated to the quantity AS*. 

To be more specific, let us write {ASy in terms of ((5q, 6p), 

f f 

i=l j=l 

+ {{QiPj)^ - {Qi)^{Pj)^)SPiH + i{Pi^j)i, - {Pi)^{i3)^)^^i^Pi ' (10) 

or 

{ASf = Sp^-f^^p + Sq^Y^Sq + 5p^7''^5q + Sq^Y'^P , (H) 

where we have introduced the matrices 7^'^ = (AiBj)^ — {Ai)^{Bj),^, and a"^ denotes the 
transposed of the vector a. To gain some insight into the meaning of this quantity, we will 
consider (AS)"^ for the one-dimensional case, 

(AS)^ = {agSpf + {ap5qf + {{qp + pq}^ - 2{q)^{P)^)SqSp , (12) 

where cXg and Cp are the root-mean-square deviations of q and p respectively. To continue 
with our discussion, a rotation in phase space is made, so that the term {{qp + pq)^ — 
2{q)^{p)^) in the previous equation is zero, and A^ is given, in terms of the new phase 
space variables, by 

{ASf = {a,5pf + (Mg)^ (13) 

where cTg (cTp) gives the support of the state in the variable q {p). A classical action A = a^ap 
can be associated to the state of the system. Eq. (fT^ implies that displacements such that 
aq6p ^ h or ap6q ^ h give AS* ^ h, and the main point of our analysis is that they also lead 
in general to a significant variation of |C^P, irrespective of the value of action A. In this 
sense values of the order or larger than h of the AS'-action scale are always needed for this 



6 



environmental system to induce decoherence. It is possible to define other relevant quantities 
with units of action. For instance, values of AZ = 6qSp leading to a significant decrease of 
the overlap are related to the size of the structure of the distribution associated to the state 
in some particular phase space representations 0]. For the displacements discussed above 
AZ ^ fi-'^/A, and if A » h the result that sub-Planck displacements on the AZ-action 
scale are relevant for the decoherence process induced by S comes naturally. Coming back 
to the multi-dimensional case, when the dimension of the problem increases more terms will 
contribute to {ASy in Eq. (fTUj). and smaller displacements in each variable are needed to 
reach the threshold AS ~ h, leading to the result that a larger number of degrees of freedom 
will favour the decoherence process 

To illustrate the difference between AZ- and AS-action scales we consider a general 
Gaussian state in one dimension 

^(g) = ( ^] ' e^Po9Ae-(9-5o)V-^ (14) 
\t{\zY ) 

with z = z^^^ izi^ zn = {h/apY , and zi = zr-^J Aa^a^ — h^/h. Straightforward calculations 
lead to the exact expression 

\C46q,6p)\' = exp[~iAS)yn'] , (15) 

where 

{ASr = {<yM? + i<^M + 1^ ^^^%^'''' )' - iSpSq)\ (16) 

in terms of the first two moments of 5*, as expected for a Gaussian wavefunction. Eq. (fT3)) 
shows that values of AS ^ h are needed to obtain a significant decrease of the overlap 
|C^p. If we now choose, for instance, particular values of the widths ag and dp so that the 
Gaussian state is much narrower in coordinate than in momentum space, say ag 
and CTp ^ (in arbitrary units), it is clear that a displacement {6q, 5p) = {y/h/2, v^/2) 

will take the shifted Gaussian completely away from the initial one. The different actions 
associated to that same displacement are AS ~ 5h and AZ = h/A, corresponding to over- 
Planck and sub-Planck values respectively. 
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III. SUB-PLANCK STRUCTURES IN PHASE SPACE DISTRIBUTIONS. 



The behaviour of the overlap |C^P with ((5q, Sp) can be alternatively studied through the 
distribution associated to the state in different phase space representations. In this section we 
will derive the relation between the overlap and the action AS* using a wide class of quantum 
quasi-probability distributions F{q, p;x) 3) the Wigner |^ and Husimi l3] functions being 
nothing but particular cases. The choice among the F functions associated to the same 
quantum state of a system, or, equivalently, the selection of a particular representation 
(given by function x), is similar to the choice of a convenient set of coordinates 12, 13, 1^ . 
Within this framework, the expectation value of any operator G{c[, p) is written as the phase 
space integral 

(G'(q,P))^ = J d^qd^P F^ici,P;x) 9 i(l,P;x) , (17) 
where F^(q, p; x) is obtained from the quantum state \ip) as 

Th 



F^{^,p; x) = ^^^yf j dfedWdfux{0,T)(u.+ '^ 



u 



-i[0-(q-u)+x-p] 



18) 



The Wigner and Husimi functions, for instance, are obtained by replacing x(0, r) = 1 and 
x{0, r) = exp{ — ^[(tA)^ + (0/A)^]} respectively. The function qfq, p; x) is the image of the 



operator G in phase space according to the kernel function x 



^(q,p;x) 



—Y [ d^ed^rdk—^, 
2vr; J x{0,t) 



u 



Th 

~2 



G 



Tfi 



j[6»-(q-u)+x-p] 



(19) 



and it is not necessarily equal to the classical magnitude. In particular, the expectation 
value of the displacement operator D can be written as the phase space average 



C^(5q, 5p) 



d^qdfp F^(q, p; x) c?(q, P; x) 
1 



d^qdfp F^{q,p;x) e*('l■'^P+P•^'^)/^ 



(20) 



x{,5p/h,5i\/Ti) 

where the function (i(q, p; x) was obtained by integrating the r.h.s. of Eq. ()19p with G 
replaced by t). Notice that Eq. (jSj) is a particular case of Eq. (j^UI) for which the Wigner 
function has been chosen as the distribution associated to the state, C^(5q, 5p) being equal 
to the Fourier transform of W^. Eq. ()2()|1 can be used to understand the relation between 
and AS* from the point of view of the phase space distribution F^. On one hand, if 
the exponential factor does not vary significantly over the support of -F^(q, p;x)) which 
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occurs for small enough values of (5q and 5p, then the overlap will only differ slightly from 
the normalisation integral of the original distribution, J d^qd^p F^(q, p;x) = 1, leading in 
general to a small decrease of the function |CJf. (Notice that x(0,0) = 1 is needed to 
guarantee that -F,/,(q, p; x) is normalised to one 13[.) The condition for these variations not 
to be significant is equivalent to the condition that the root-mean-square deviation of S, 
AS, is smaller than h. On the other hand, to obtain significant decay of the overlap, rapid 
oscillations, with AS at least of the order of h, are needed. Due to the properties of the 
Fourier transform, and since the value of xi^p/^, ^q/^) is close to one for small enough 5q 
and 6p, the initial decay of with the displacement is related to the large scale structure 
of the distribution F^, that depends mainly on the size of the state support in phase space. 
However, the detailed behaviour of |C^P with arbitrary displacements, and, in particular, 
qualitative features like oscillations, will depend on the state under study. 

A different question is how sub-Plank structures emerge in some phase space distributions 
associated to the state and how they are related to the main features of Cw,. For kernel 

n 

functions such that \x{6,t)\ = 1, the corresponding distributions verify [1^ 

|C^(<5q,5p)|' = {2TThy J dfqdfpF4q,p)F4q + 5q,p + 5p). (21) 

(Notice that Eq. (jD) is a particular case of Eq. (jH)).) For these representations, the fact 
that a given displacement leads to |C^P ^ is manifested in a complex structure of the 
distribution F^ on the scale of the AZ-action for that displacement, as pointed out in Ref. 
P] for the Wigner distribution. When displacements with AZ << h lead to small values of 
|C^P, the distribution F^ will show a complex structure at sub-Planck scales. This result is 
a consequence of the particular choice of the phase space distribution. For the same state, 
the Husimi distribution (obtained by smoothing the Wigner function, so eliminating the 
sub-Planck scale structure) will lead to the same overlap |C^p. This is not surprising since 
the overlap depends on the state, and that dependence is manifested in different ways for 
different phase space representations. 

IV. A ONE-DIMENSIONAL TIME-DEPENDENT ENVIRONMENTAL SYSTEM 

In this section we will analyse the dependence of the overlap |C^P on the action AS{6q, 6p) 
in the context of a particular one- dimensional model for the environmental system S, de- 



9 



scribed by Hamiltonian 

= K COS (g — I sin t) + -^o.(f ■ (22) 

This quantum model system has been previously used in the context of decoherence 
and describes a particle of mass m = 1 (arbitrary units are used throughout) confined by a 
harmonic potential that is perturbed by a spatially and temporally periodic term. For the 
parameter values used in this work, k, = 0.36, a = 0.01, / = 3.8, and h = 0.16, the motion 
in the classical counterpart of this system exhibits a chaotic character 

To prepare the state of the environmental system prior to the interaction, we let a given 
initial state evolve until preparation time T, when it is coupled to the pointer system S. 
The coupling strength is assumed to be high enough so that the two-systems evolution can 
be followed as described in the introduction, i.e., by neglecting any contribution coming 
from the dynamics induced by Hamiltonian (f^^ during the interaction time. This approach 
allows us to discuss the values of actions involved in the decoherence process induced on 
system S in terms of general displacements in phase space, irrespective of the detailed values 
of the coupling constants and interaction times [isj ]. 

As the initial state {T = t = 0) for the preparation process we have chosen a coherent 
state la) = e"'"'^''^ J2'^=o 7^ l'^) harmonic oscillator Hqa = p^/(2m) +ag^/2, where 

\n) is the eigenstate of Hqa with energy {n + l/2)fi^ya/m. (We have checked other possible 
initial states obtaining qualitatively similar results.) Its time pro pag ation under Hamiltonian 
(I22|) has been obtained by means of the split-operator method 19|. 

As time T increases, the state spreads in coordinate as well as in momentum space through 
the available phase space as shown in the insets of Fig. To characterise this dynamics 
the quantity 

a = ^ (23) 

is used (see Fig. P). It shows a rapid initial decay (until time T ^ 20), followed by a much 
slower decrease for longer times. The behaviour of a for small preparation time T is related 
to the fast initial increase of the widths and Up. The variation of a for longer times is 
mainly due to the time-dependent term in the Hamiltonian. Should not be for the presence 
of this time-dependent term a would not decrease beyond a certain minimum value related to 
the maximum position and momentum widths compatible with a fixed mean system energy. 
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FIG. 1: Action a = fi^ /{aqap) as a function of the preparation time T (solid line) for the initial 
coherent state with a = f>i. The action AZq needed for a displacement in the direction 5q ~ Q.SSp 
to reduce the value of |C^P to 0.5 is also shown (dashed line). The inset shows the dependence of 
the widths Cg (left) and Cp (right) with the preparation time. (Notice the different scales in the 
vertical axis for each case.) Arbitrary units are used. 

Fig. 121 shows |C^P versus AS* for different preparation times and for a given direction 
in phase space. (The results for any other direction show the same qualitative features.) 
The different curves, corresponding to different preparation times, decay in the same AS- 
scale. To emphasise this result we represent in the inset the value of AS* needed to obtain 
|C^P = 0.5 versus T. The AS-action values for any preparation time are of the order of 
supporting AS* ^ as a relevant scale for the studied decoherence process. (Notice that the 
apparent convergence of A So to a value close to % is only a consequence of the chosen value 
for \C^\\) 

The values of AZ for which |C^P = 0.5 are also shown in Fig. G]and Fig. |2l After some 
time T ^ 5, the action a sets the scale of the random structure developed in the distribution 
associated to the states of this system in some phase space representations, for example in 
the Wigner function Taking into account the discussion below Eq. ()21|). the action AZ 
for displacements producing a significant decrease of |C^P will be of the order of the action 



after T 5, as shown in Fig. ^ 
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FIG. 2: |C^P as a function of the action AS{6q, 5p) in the direction Sq ~ 6.86p in phase space 
and for different preparation times: T = (sohd hne), T = 10 (dashed hne), T = 20 (dotted hne), 
and T = 500 (circles). (Same initial state as in Fig. The inset shows the actions ASq (upper 
curve) and AZq (lower curve) needed for the displacement to reduce |C^P to the value 0.5 versus 
the preparation time T. The straight line shows the action h for reference. Arbitrary units are 
used. 

V. NON-LINEAR CONFINED ENVIRONMENTAL SYSTEMS 

For now on, a different model for the environmental system S will be considered, that of 
a time-independent Hamiltonian, Hnl, with a non-linear confining potential and a discrete 
energy spectrum. Under certain assumptions, this model will allow us to obtain analytical 
expressions for the overlap between the states 14^+) and {ip-)- Instead of considering a 
particular environmental state obtained after some fixed preparation time T, we will 
study the dependence of the overlap averaged over the preparation time on an averaged 
AS-action. For non-linear confined systems, the main features of this stationary description 
can be associated to all the states prepared from a given initial one, provided that their 
preparation time is long enough. In the first part of this section we will determine the 
stationary properties of relevant to our discussion. We will assume that the states are 
prepared from a given \ip{T = 0)), and make use of the Wigner distribution in phase space 



12 



associated to them. Although the procedure and the results are independent of the choice 
of a particular phase space representation, the use of the Wigner distribution will allows us 
to extend our analysis afterwards for systems for which the Berry- Voros conjecture is valid. 



A. Stationary properties of the overlap 

In the basis of eigenstates of the Hamiltonian Hnl-, which will be assumed to have, for 
simplicity, a non-degenerate spectrum, the wave function at preparation time T is given by 

V^(q,T) = ^c„e-^^"^/V(q), (24) 

n 

where H(fn{^) = En(fn{^) and c„ = / (ig (q)'?/'(q, 0). The Wigner distribution is obtained 
introducing expansion into Eq. Q. Splitting the result into time-independent and 
time-dependent terms, 

w4ci,p,T) = 5^|c„|Xn(q,p) (25) 

n 

where W^^{q,p) is the Wigner distribution associated to the energy eigenstate v^n(q)- For 
non-linear systems, it turns out that the Wigner distribution spreads from its initial (T = 0) 
support in phase space until it occupies most of the available phase space volume at some 
preparation time Tc. From time Tc on, the small details of the Wigner distribution will 
change with time, but in general its long scale structure will remain as a stationary property. 
To extract that characteristic long scale structure we employ the time-averaged Wigner 
distribution 



;q,p) = lim- [ dTW4q,p,T) = Y^Kl^W^M^P) , (26) 



where we have taken 
Introducing into Eq. (0), we obtain the time-averaged quantity 



^limi^'rfr^c„c>-^(^"-^-)^/^| ^^e^'''P/Vn(q-qV2)<^:;.(q + q72) = 0. (27) 



C^{Sq,5p) = j d^grf^pe^(P-^''+^-^P)/"Wv(q.P). (28) 
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that describes the stationary properties of the overlap between \'4)+) and \ip-)- According to 
Eq. (PHj) . can be identified as the generating function of all moments of 5 = p • 5q+q ■ 5p 
with respect to the distribution W^. Therefore a set of equations similar to Eq. and 
(jHI can be obtained. These equations imply that the initial decay of is ruled by the 
fluctuation properties of S at stationary conditions. The action scale involved in the decay of 
for small displacements (5q, (5p) can be in general associated to any state with preparation 
time longer than T^. This result follows from describing properly the long scale structure 
for T > Tc and the discussion in Sec. III. In the rest of the section we consider a family of 
quantum systems for which can be obtained analytically. 



B. Systems described by the Berry- Voros conjecture 



We shall now pay special attention to (1) quantum systems with time-independent Hamil- 
tonians and classical chaotic counterpart and (2) regular quantum systems with particular 



random components in their potentials 



2ll |. for which the relevant quantities are obtained 



after averaging over the noise. There are both, experimental and numerical evidences, that 
for these systems the so called Berry- Voros conjecture is valid, namely, that one can approx- 
imate the Wigner density associated to an energy eigenstate by a microcanonical density 



23 



24| 



1 6{E^-H{ci,p)) 



(29) 



/ (2^^)/ ^ ~ -^(ct, p)) is the local average density of states at enep 



2G, m 



{2nK)S p{Er, 

where p(E„) = / f^' ...... 

En, and H{q^, p) is the classical Hamiltonian associated to the quantum one j9|, 
l2fil |. The Wigner distribution in the semiclassical limit fills the available phase space that 
corresponds to an energy shell of thickness of the order h and its amplitude fluctuates around 
the microcanonical density. Furthermore the density function in Eq. (j29|) is just the leading 
approximation of a semiclassical expression for W^(^„(q, p). The next to the leading terms 
depend on the periodic orbits of the classical system and take into account the possible scars 
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28 



3- 



Replacing M^(p„(q, p), implicit in Eq. (j28|) . by the expression in Eq. (j29|l . it follows 

' \ / iPn 



(30) 
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where 

' e 



is the microcanonical average of e*"^/^. For a Hamihonian of the form if(q, p) = p^/2M + 
\^(q), and after integrating over the momentum variables, one obtains 

xJi_i(^V2M(E„-y(q))), (32) 

where is the Bessel function of order //2 — 1. Eqs. ()3()|1 and ()32|1 lead to a formal 

^ -By 



expression of the time- averaged two-point correlation function C^(5q, 5p) in terms of 
the potential V{q). These equations constitute the main result of this section and are the 
starting point for the analysis of particular examples. In the following we shall particularise 
Eq. ()30p for systems with a random component in the potential such that the average over 
the noise is the /-dimensional harmonic potential. 



1. The f -dimensional harmonic oscillator 



The classical Hamiltonian for a generic /-dimensional harmonic oscillator, 

2m,- ■ 2 



can be rewritten, in terms of the rescaled coordinates and momenta 

fW 
Pi = \ — Pi 

V rrii 



(34) 



as the spherical harmonic oscillator 



^(q, P) = 7^ipl + Pl + ---+p'f) + ^Mu'^iql + ql + ... + qj). (35) 



2M 

The integral in Eq. (jHTj) reads 
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with 



After some manipulations, it follows that 



and integrating over variable ^, 



BV . ,,,/ 



'i(p-5q+q-(5p)/n\ _ 2/ 



(39) 

For an eigenstate of energy p{En) = El~^ / {T{f){hujy) , where r(/) denotes the Gamma 
function of argument /. Besides, „ = (p^) = ME„, a^^ = (q^) = En/Mu^ (for this 
case the mean values of position and of momentum vanish), giving 



^{p.5q+q.5p)/n\ ^ 2^-1 T{f)p{E, 



-p-|- f I \^ / r \ / j 

iU=l^i/ /^ l'?Pkg,n -|2 I j- |^qkp,n -|2yf-l 

where the characteristic action 

A^f = ^J{\5v\<y,,nf + {\H<yv,n? (41) 

has been introduced. AS*^^ is nothing but the action AS* introduced in Eq. (fTUj) calculated 
for the nth eigenstate using the Berry- Voros conjecture. For the superposition state (|^. 
one obtains 
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FIG. 3: p versus AS" (dashed line) for the direction 6q ~ 6.86p in phase space and the 

same initial state as in Fig. ^ For comparison, |C^P is shown for the initial (solid line) and 
T = 500 (dots) states of Fig. |21 Arbitrary units are used. 

so that the typical action that controls the decay of the overlap is the one related with the 
coefficients c„ that contribute more to the initial state. The action AS" , evaluated for the 
average distribution under the Berry- Voros conjecture, is related with (AS"^^) by 

n 

For any displacement (5q, 5p) the previous relation can be inverted and used to write 
in terms of AS*^^. 

To illustrate this result, in Fig. 01 we plot versus AS* for the one-dimensional case 

/ = 1. The AS'-action scale for the decay of the overlap is dictated by the value h = 0.16, 
as in the case described in the previous section. Note that this result is expected as a power 
expansion of the Bessel function Jq in Eq. ()42j) to the second order in AS" consistently 
recovers the result in Eq. (jHJ. 
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VI. DISCUSSION 



The results of previous sections show that the relevant AS'-action scale to the decay of the 
overlap |C^P for small displacements (5q, 5p) is given by h. The one- dimensional Gaussian 
state is a special example for which the dependence of on AS* is given explicitly by 

Eq. ()15|) . and its monotonic exponential decay is independent of particular details of the 
state, as for instance the widths in position and momentum. (On the contrary, the decay of 
the overlap with the displacement will depend on (jq and dp through Eq. (fT^ V In Figures 
121 and El the exponential dependence associated to an initial (Gaussian) coherent state is 
compared to the one corresponding to states at different preparation times. Although all the 
curves shows a similar initial decay, (dictated by AS* ~ h), the ulterior behaviour can have 
qualitatively different features, the presence of oscillations in the overlap for intermediate 
values of AS" being the most relevant one. It is worth noting that these oscillations can 
never be regarded as true revivals. |C^P can be interpreted as the overlap between the 
states \ip) and D((5q, 5p)| ■?/;), the second one being obtained by a rigid displacement of 
This implies that |C^P can not be equal to one for non-zero displacements since the support 
of the state in phase-space in finite. However, large amplitude oscillations are possible as 
shown in Fig. |21for T = 10. 

The pattern of oscillations will change in general with the preparation time. In the system 
described in Sec. IIVI no oscillations are present for the initial state. For small preparation 
times some oscillations appear (see Fig. |21for T = 10) but their amplitude decrease when 
the preparation time increases. For larger preparation times only oscillations with small 
amplitude are found. This behaviour can be interpreted by using Eq. (j21|) with the Wigner 
function, as proposed in Ref. j^. For T = 0, the Wigner distribution associated to the 
coherent initial state is a Gaussian in phase space, and the monotonic decrease of |C^P with 
AS reflects the decrease of the overlapping regions between the states \ip+) and lip^) (or, 
equivalently, between {tp) and D\iIj)). For small preparation times, the isolated evolution of 
the environmental system prior to the coupling generates a regular large scale structure in 
the distribution (characterised by large values of AZq in Fig|2I). For this case the coincidence 
between maxima and minima of that large scale structure in and {ip-) is responsible 
for the oscillations in the overlap. For longer preparation times, smaller scale structures 
appear in the distribution (corresponding to smaller values of AZq), and more importantly, 
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the randomness of the distribution of the patches in the structure increases (reflected in the 
similarity of the actions AZq and a). Then, as T increases the amplitude of the oscillations 
becomes smaller until they are eventually negligible. This behaviour is expected in general 
for any non-linear system, with the only difference in the preparation time T needed to 
develop the small scale structure. 

In the light of this discussion, special care must be taken in the interpretation of the results 
of Sec. where broad oscillations in the time averaged overlap \C^ p could appear for 

^BV 



large AS" (see Fig. As the systems considered are non- linear, the states will develop 
in general a complex small random structure for long enough preparation times, and only 
negligible oscillations will be present in the overlap. The broad oscillations in \C^ p 
are the result of the use of the Berry- Voros conjecture, that describes correctly the large 
scale structure but fails in describing the small scale correlations. Therefore, following the 
discussion in Sec. IIIII related to Eq. (pUj) . only the initial decay (corresponding to small 
displacements) for each particular sufficiently long preparation time is well described by 
l^lp I ■ 

In the approach used in this work, the effect of the coupled evolution in the environ- 
mental system is equivalent to rigid displacements in phase-space of the state \iIj{T)) to 
give \ip+(T;6t)) and \ip-(T;6t)). (The dependence of \4'±) with the interaction time 6t is 
made explicit.) No additional structure in phase-space in the states \'ip±) is generated during 
the coupling, as the contribution of Hs is neglected. The interaction time 6to required to 
obtain a value |CoP of the overlap is given by the condition |C^(5qo, (5po)P = jC'oP, where 
the magnitude of the displacements are 5qo = — 2Cp5to and 5po = — 2Cq(5to- Therefore the 
larger the coupling constants, the smaller the interaction time SIq. The condition AS* ~ h 
establishes a lower bound for the value of the displacements and consequently for the inter- 
action time needed to attain effective decoherence. An alternative derivation of the lower 
bound is pointed out in Ref. A different aspect is the dependence of this 6tQ with the 
environmental state prior to the interaction. The size of the displacements (5qo, 5po) can 
be described by the action AZq = 6qo6po for each particular state. As discussed in Sec. 
mil AZq is of the order of the action that sets the scale of the structures in the distribution 
for some phase space representations fulfilling Eq. (PT|) . As a result, 6tQ decreases as the 
structure in the distribution associated to the state becomes smaller. For example, in the 
system analysed in FiglH the interaction time 6to is proportional to i/a provided the 
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preparation time is long enough for a to describe properly the small scale structure. A more 
complex situation appears when the evolution induced by is not neglected 

0Q- In 

that case, besides the displacement, the distribution of the structure in phase space of the 
states and will change during the interaction. For the system studied in FigQ 
two different regimes can be distinguished. For T ^ 20, a rapid variation of the sizes of the 
structure with time is found and both mechanisms, the displacement and the development of 
structure, will determine the interaction time 5to- However, for T ^ 20, the variation of the 
sizes of the structure is much slower and dt^ is determined by the time required to produce 
the displacement in phase space. As the displacement is approximately independent of the 
details of the state, 51^ will be weakly dependent on the preparation time for T > 20 

Another important point to discuss is the dependence of the decoherence process with 
the number of degrees of freedom of the environmental system. As the number of degrees of 
freedom increases, smaller displacements in each variable are needed to obtain AS" ~ that 
sets the action scale for the initial decay of the overlap in all cases, and the corresponding 
interaction time will be smaller too. This is compatible with the observation that the larger 
the environment the more effective the decoherence process. 

Experimental tests of the decoherence process in the context discussed in this work can 
be in principle realized in the systems described in Refs. ||3(J]. The interaction between two 
oscillators is mediated by a term of the form hG aj^ag (a^+a^), corresponding to a scattering 
process in which a quantum of energy of the environmental system S can be absorbed (a^) 
or emitted (a^.) whereas the number of quanta of the pointer system S remains the same. 
For these cases, the coherences of the reduced density operator of the pointer system in the 
basis given by the Fock states are proportional to the overlap between the states D{a = 
iGn6t)\ilj{T)) and D{a' = iGn'5t)\ilj(T)) . The operator D[a = iGndt) = expjaa^ — a*as} 
produces a displacement in phase space that depends linearly on the interaction time 5t, 
the coupling constant G, and the index n of one of the Fock states of S involved in the 
coherence under consideration. 

In summary, the role of as a boundary between different decoherence regimes has been 
clarified in the context of a characteristic action AS, which depends on the quantum state 
of the environmental system. We related the action AS* with the complementary quantity 
AZ, and described their connection with the pattern of structures developed in phase space. 
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